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Recently we have used spectroscopic mapping with the scanning tunneling microscope to probe mod- 
ulations of the electronic density of states in single crystals of the high temperature superconductor 
Bi2Sr2CaCu20 8+ 5 (Bi-2212) as a function of temperature [C. V. Parker et al., Nature (London) 468, 
677 (2010)]. These measurements showed Cu-0 bond-oriented modulations that form below the pseu- 
dogap temperature with a temperature-dependent energy dispersion displaying different behaviors in the 
superconducting and pseudogap states. Here we demonstrate that quasiparticle scattering off impurities 
does not capture the experimentally observed energy- and temperature-dependence of these modulations. 
Instead, a model of scattering of quasiparticles from short-range stripe order, with periodicity near four 
lattice constants (4a), reproduces the experimentally observed energy dispersion of the bond-oriented 
modulations and its temperature dependence across the superconducting critical temperature, T c . The 
present study confirms the existence of short-range stripe order in Bi-2212. 

PACS numbers: 74.72.Gh, 74.55,+v, 74.20.Pq 



Understanding the interplay between superconductivity, 
the pseudogap [lj, and the possibility of charge and spin 
ordering phenomena in the proximity to the Mott insulat- 
ing ground state continues to be one of the most challeng- 
ing problems in condensed matter physics. Periodic pat- 
terns of spin and charge, referred to as stripe order [EH6), 
have been detected in static form in the La-based cuprates 
0. Yet their relevance to the wider class of high-T c com- 
pounds and significance to the mechanism of superconduc- 
tivity and pseudogap remains unknown. Addressing this 
question is challenging since in the presence of disorder, or 
in the absence of favorable structural distortion (71 El, such 
an order is likely short range or fluctuating and hence hard 
to detect l9lfT0l. While the observation of static stripe pat- 
terns with scattering techniques is well established, detect- 
ing fluctuating order and distinguishing it from other elec- 
tronic spatial modulations is still being developed [9]. Es- 
tablishing the presence of fluctuating or short-range stripe 
order in the cuprates and understanding its correlation with 
other phenomena, such as superconductivity or the pseudo- 
gap, is of great importance. 

In this paper we analyze spectroscopic measurements 
with the STM that allow us to probe the spatial variation 
of the local density of states (LDOS) in single crystals of 
Bi-2212. Recently, we have shown that spatial features of 
the LDOS associated with incipient stripe order ifTTTl can 
be distinguished from those due to impurity-induced quasi- 
particle interference H2HI3 in this compound. These pre- 
vious measurements demonstrate that signatures of incip- 
ient stripe order, which appear as bond-oriented modula- 
tions in the LDOS |[T6l - [T9ll , first appear below the pseu- 
dogap temperature T* across a wide range of doping. We 
have also shown these modulations in Bi-2212 to have the 



strongest intensity when the samples hole doping is close 
to 1/8, the concentration at which static stripes have been 
found in La-based cuprates [3 1 and as predicted by most 
stripe models [4-6, 9]. Yet, the wavelength associated with 
these modulations shows a significant energy-dependence, 
a behavior not expected for static long-range order. Fur- 
thermore, the energy-momentum structure of these modu- 
lations in near optimally doped samples exhibits dramati- 
cally different behaviors across the superconducting transi- 
tion temperature T c . Here we show that a model of scatter- 
ing of quasiparticles from short range stripe order can cap- 
ture the energy-dependence of the LDOS modulations both 
above and below T c . This model also provides an insight 
into the particle-hole symmetry of the electronic states in 
the superconducting and pseudogap phases. These findings 
demonstrate the strong interplay between incipient stripe 
order and the electronic properties of both superconducting 
and pseudogap states of the cuprates. Establishing this con- 
nection in compounds other than the La-based compounds 
suggests that incipient stripe order plays an important role 
in all families of cuprates. 

Figure [TJa) shows an example of real space mapping of 
the low-energy conductance (G(r, V) = dI/dV(r,V)) 
using STM on an optimally doped (OP91, T c = 91 K) 
Bi-2212 sample, carried out at 35 K, below T c . Discrete 
Fourier transforms (DFT) of such conductance maps show 
strong peaks at wavevectors marked Q* along the Cu-0 
bond direction (Fig. |TJb)). As we have previously re- 
ported, the Q* peaks are present in the DFT maps above 
T c and up to the pseudogap temperature T* IfTTTl . De- 
spite the lack of sensitivity of the Q* intensity to the on- 
set of superconductivity IfTTTl . we find these modulation to 
show systematic changes in their energy dependence with 
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FIG. 1. Optimally doped Bi-2212 sample (T c = 91 K). (a) Real space mapping of the low-energy conductance, G(r, V = 10 mV), 
taken at 35 K, and measured at a setup bias of —150 mV and setup current of 40 pA. (b) DFT of the data in (a), showing strong peaks 
corresponding to the Cu-0 bond oriented modulations marked Q* (red arrow) and b-axis supermodulation (green arrow), (c) DFT 
of Z(r,V = 10 mV) on the same area and temperature. Inset shows schematic of the octet model predicting seven wavevectors (qi 
through qj), four of which can be seen in the DFT in (c). Corner of the DFT maps (a,b) correspond to (27r, 2tv) in units of 1/ao, where 
ao = V2a and a is the nearest neighbor Cu-Cu distance. PSD stands for power spectral density. 



temperature. As shown in Fig. |2ja,b), for both optimally 
doped (OP91) and underdoped (UD84, T c = 84 K) sam- 
ples, the Q* wavevector disperses with energy, showing 
a particle-hole symmetric dispersion at temperatures well 
below T c , with the Q* wavevector being approximately 
symmetric between positive and negative sample bias. At 
higher temperatures, specifically for T > T c , we see that 
the energy dispersion, Q*(E), is no longer particle-hole 
symmetric. While the Q* modulation is present through 
the entire temperature range below the pseudogap temper- 
ature T*, the particle-hole symmetry (asymmetry) of its 
dispersion correlates with the presence (absence) of super- 
conductivity below (above) T c . 

In the superconducting state, well below T c , we 
find that ratios of conductance maps Z(r,V) = 
G(r,+V)/G(r,—V) fT4l provide an effective way of 
enhancing features of the data associated with supercon- 
ductivity. As previously demonstrated these features of 
Z-map (marked as g's, Fig. [TJc)) are similar to those 
expected from the so-called octet model for scattering 
of Bogoliubov-de Gennes quasiparticles (BdG-QPI) from 
random impurities in a d-wave superconductor (schemat- 
ically shown in inset of Fig. QJc)) EMS El EE For 
temperatures well below T c , Q*(E) from the Z-maps ex- 
hibits energy-dispersion trends that are in poor quantita- 
tive agreement with the the q 1 vector of the octet model 
extracted from the band structure obtained from angle- 
resolved photoemission (ARPES) measurements (9l [22) 
(dashed gray line in Fig. |2jc,d)). However, by increas- 
ing the temperature above T c , where Q*(E) is no longer 
particle-hole symmetric in the conductance maps, we find 
an artificial flattening of Q* (E) in the Z-map analysis (Fig. 
[2jc,d)). This is evident since use of the Z-map assumes 
that the tunneling excitations measured by the STM tip are 
particle-hole symmetric BdG quasiparticles [23], a condi- 
tion true only in the superconducting state (Fig. |2ja,b)). 
Fluctuating superconductivity also cannot explain the be- 
havior of Q*{E) above T c , since such scenario would 
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FIG. 2. Q*(E) for OP91 samples (a) and UD84 samples (b) 
at selected temperatures, measured along the Cu-0 direction in 
units of 27r/a, where a is the nearest neighbor Cu-Cu distance. 
Q*(E) extracted from Z-map analysis for OP91 samples (c) and 
UD84 samples (d), for the same temperatures as in (a) and (b) 
respectively. Octet model prediction for q\ is displayed as a gray 
dashed line in (c,d). The experimental data are reproduced from 
results reported in |[TT1 . 



as well require the presence of particle-hole symmetry in 
Q*{E) in clear disagreement with our observations (Fig. 
[2]). Therefore, to determine the nature of the Q* modula- 
tions not only below T c (particle-hole symmetric Q*(E)) 
but also above (particle-hole asymmetric Q* (£")), we fo- 
cus only on information obtained from the conductance 
maps and their temperature dependence for the rest of this 
paper. 

To obtain a unified understanding of the Q* modulations 
as a function of temperature, we calculate the quasiparticle 
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scattering from a single impurity with a T-matrix approach 
EUlh in both the superconducting and the normal states, 
and consider how the results of such scattering would be 
modified in the presence of short range stripe order Il9ll24ll. 
In this approach, the LDOS modulations due to impurity 
scattering are initially determined by 

Xo(q^) = ■^-Im[A 11 (q,(j) + A 22 (q,-u)] (1) 

where A(q,u) is a 2 x 2 matrix, as prescribed by the 
Nambu-Gor'kov spinor formalism: 

d 2 k 



(2tt)2 

For a single impurity, the T-matrix is 



Go(k + q,u)T(u)G (k,u) (2) 



T~\uj) = (V s a 3 + Vml)- 1 - J ^G (fc, 



u) (3) 



where G is the unperturbed (2 x 2) Nambu-Gor'kov single 
particle Green's function and V m (V s ) is the spin-flip (non- 
spin-flip) component of the impurity potential. 

We assume the simplest form for the Green's function, 



G 1 (fc, u) = (cu + iT)I - e k a 3 - A fe ai 



(4) 



where e k is the band structure obtained from photoemis- 
sion experiments [22 j, T (3 mV) represents quasiparticle 
broadening, I is the identity matrix, G\ and cr 3 are Pauli 
matrices, and A k is the d-wave superconducting gap func- 
tion. Above T c , where A k = 0, we find that the results of 
the above calculations, which are insensitive to the choice 
of the impurity potential (V m , V 8 ), do not produce any 
peaks in the vicinity of Q* wavevectors within the experi- 
mental energy range (-40 mV to +40 mV) (see for example 
Fig. |3ja)). Instead, the square-like contours seen in Fig. 
[3ja) persist for all energies, although with varying inten- 
sities and wavelengths, in clear disagreement with the data 
(see for example Fig. [3Jc)). This behavior indicates that the 
modulations near Q* are not solely due to impurity scatter- 
ing, rather as we have previously argued, a consequence of 
short-range stripe order in Bi-2212. 

To model the scattering from incipient short-range stripe 
order, we follow an approach proposed by Kivelson et al. 
0. The presence of a short-range stripe order centered 
at Q st , corresponding to a periodicity of 4a, strongly en- 
hances the quasiparticle scattering with momentum trans- 
fer near Q st . While impurity scattering alone (% ) might 
not produce peaks near Q st , scattering from stripes will 
enhance any weak momentum- and energy-dependent fea- 
tures of Xo n ^ar Q st . Therefore, the LDOS measured by 
STM (Xdos{Qi ^0) will exhibit an energy dispersing peak, 
as long as Q st has finite correlation length (short-range). 

The simplest phenomelogical way to model Xdos m the 
presence of scattering from this incipient stripe order is to 
write it as 
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FIG. 3. (a) xo(q,u = 10 mV) calculated for A k = 0. (b) 
Xdos(q^) calculated from Eq. fusing xo(<?,^) in (a), (c) 
DFT of G(r, V = 10 mV) on an OP91 sample at T = 35 K 
(setup bias and current equal —150 mV and 20 pA, respectively), 
(d-g) xo(Qi w) calculated for A k = 45 meV at the antinode, for 
selected energies, (h-k) xdos(q^) calculated from xo(q^) 
in (d-g) respectively, (l-o) DFTs of G(r, V) on an OP91 sam- 
ple at T — 35 K, for selected energies (setup bias and current 
equal —150 mV and 40 pA, respectively), (p) Schematic rep- 
resentation of Eq. [5] where xo(q,u = mV) is calculated in 
the superconducting case, (q) Line cuts through xo(q, w) (blue), 
Xst(q) (green), and xdos(q,v) (red) along the direction of Q* 
(see lines in (p)). To obtain figures (b, h-k, p) Xst(q) was calcu- 
lated from Eq. [S]for f = 0.1(2?r/a), and R = 3. (r) Similar line 
cuts for uj = 30 mV (xo displayed in (g) and xdos m (k)). 



where the stripe induced density of states takes the form 

(9 - Qst) 



Xst(q) = Rx exp 



(6) 



Xdos(q^) = Xo(q) x (Xst(q^) + 1) (5) where Q st = (±0.25, ±0.25)27r/a , £ determines the in- 
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FIG. 4. Comparison of Q*(E), above and below T c , obtained 
by the same fitting procedure for the model calculation and the 
experimental data measured for OP91. Q*(E) is measured in 
units of 27r/a where a is the nearest neighbor Cu-Cu distance. 

verse correlation length and R is the scattering enhance- 
ment at Q st . This choice of ordering wavevector is also 
supported by the doping dependence of the Q* intensity we 
have previously reported ifTTIl . The result of this multipli- 
cation is a stripe-induced enhancement of the quasiparticle 
scattering near Q st (e.g., Fig. [3jb)) in good agreement with 
the data (Fig. [3jc)). Note that the strong peak in the cen- 
ter of the data is due to long- wavelength inhomegeneities, 
which are absent in the model calculation. 

The superconducting case below T c , with A k = 
A [cos(k x ) — cos(fc 2/ )]/2, yields similar results, where 
in the absence of stripe order, Xo does not produce peaks 
near Q* for low energies (—20 mV< E < 10 mV, see for 
example Figs. [3Je,f)) in disagreement with the data (Figs. 
[3jl-o)). However, in the presence of stripe order Xdos re- 
produces the peaks near Q* for all energies (e.g., Figs. [5Jh- 
k)), in good agreement with the experimental data (Figs. 
[3jl-o)). Note that the momentum structure of Xo is encom- 
passed by Xdos due to the short-range (wide) nature of 
Q st (Figs.gq,r)). 

To quantitatively extract the energy dispersion of Q* (E) 
from Xdos, we use the same fitting procedure used on the 
experimental data. Remarkably, Q* (E) extracted from our 
model calculation disperses both above and below T c and 
shows an excellent quantitative agreement with the experi- 
mental data for near optimal doping (Fig. [4]). Inclusion of 
both the dispersive feature of the hole-like band structure 
as well as an incipient stripe order appears to be required 
to capture the behavior of Q* (E). 

In contrast to the data and calculation above T c , the pres- 



ence of a superconducting d-wave gap results in particle- 
hole symmetry near zero bias below T c . Although the 
choice of scattering potential parameter does not make any 
difference in our results for T > T c , we find that below T c 
the best agreement with the experimental data is obtained 
for Vm = 100 mV, V s = mV. This choice of parameters 
is similar to the one initially used to justify the use of octet 
model below T c [20 ] and is also supported by experimental 
measurements of scattering in a magnetic field l25l . 

The success in quantitative modeling of the measure- 
ments of the energy dependence of Q* (E) modulations re- 
ported here demonstrates the importance of incipient short- 
range stripe order in understanding the electronic structure 
of the pseudogap phase and extends its relevance to the en- 
tire cuprate family, beyond the La-based compounds. In 
this model the energy dependence of Q*{E) is determined 
by the symmetries of the electronic band structure. It there- 
fore clarifies not only the presence of these modulations in 
the superconducting and normal states but also their dra- 
matic distinction across T c . The emerging physical picture 
is that although excitations of the high temperature cuprate 
Bi-2212 system have quasiparticle-like characteristics ris- 
ing from the ARPES measured band structures below and 
above T c , they also strongly scatter from a stripe-like in- 
cipient order that appears to develop below the pseudogap 
temperature T*. 
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